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The Kim-Nilles mechanism relates the (i term with the axion scale f a , leading to 
the axino-HiggsinoHiggs Yukawa coupling of order fi/f a . This can bring a dangerous 
thermal production of axinos. If the axino is stable, its mass has to be as small as 
O(O.lkeV), or the reheat temperature should be lower than O(lOGeV) taking the 
lower axion scale 10 10 GeV in order not to overdose the Universe. If the axino 
decays to a neutralino, the overproduced neutralinos can re-annihilate appropriately 
to saturate the observed dark matter density if the annihilation rate is of order 
10- 8 GeV" 2 for the axion scale larger than about 10 11 GeV. Thus, a light Higgsino- 
like lightest supersymmetric particle with a sizable bino mixture becomes a good 
c~| ' dark matter candidate whose nucleonic cross-section is larger than about 10 _45 cm 2 . 
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The strong CP problem is elegantly solved by the axion solution leading to the dynamical 
9 term: 

^ 1 n 2 n 

i a " 32tt 2 /„ ° ' 1 J 

where f a is the axion decay constant allowed in the range of 10 10 GeV < f a < 10 12 GeV 
■^j- \ In the supersymmetric extension of the Standard Model, the origin of the supersymmetric 
Higgs mass term, /i, can be related to the axion scale through the realization of the DFSZ 
axion 0]. This is the Kim-Nilles mechanism (H which introduces the /i term interaction, 

*^ ; Wkn = X^H^ (2) 

S ■ 

where Mp = 2.4 x 10 18 GeV is the reduced Planck mass. Recall that S contains the axion 
degree of freedom and generates the axion scale through its vacuum expectation value at an 
intermediate scale, (S) ~ f a , and thereby induces the \x term: \x ~ / 2 /Mp. The superpartner 
of the axion, the axino a, gets massive upon supersymmetry breaking, and can be the lightest 
supersymmetric particle (LSP) to become a dark matter candidate. However, axinos can be 
efficiently produced from the thermal bath and typically their abundance is problematically 
large. 

In the case of the KSVZ axion [H, the axino-gluino-gluon interaction, a supersymmetric 
counterpart of the axion-gluon-gluon interaction ([I]), 

^ CD = S^T a ~ aa " vrG ^ + h - c -' (3) 

is the main source for the thermal production of axinos either through freeze-out 0] or 
through regeneration [shs]. Such processes are so efficient that the axino must be very 
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light in order not to overdose the Universe unless the reheat temperature is below the 
electroweak scale to suppress its thermal production. Such a light axino requires special 
arrangements in the superpotential or supersymmetry breaking parameters to suppress tree 
and loop contributions to its mass which is generically of order the gravitino mass in gravity 
mediated supersymmetry breaking mo dels MBB. On the other hand, smaller axino mass 



may arise more naturally in gauge mediation models [Hi, 112 



So far, the cosmic axino production has been studied extensively concentrating on the 
KSVZ axion model. In this paper, we consider the supersymmetric DFSZ axion model 



realized in the Kim-Nilles mechanism extending limited discussions in Refs. [13|, [14j . The 
DFSZ axino production is led by the effective /i term interaction ([2]) which, after the field 
S gets a vacuum expectation value, induces the axino-Higgsino-Higgs Yukawa term: 

£f N = chj- a [HiH 2 + H 2 H X ] + /i.e., (4) 

J a 

where c# is a mo del- dependent parameter of order one determined by the precise relation 
between (S) and f a and the axino fraction in the fermion component S. 

As will be discussed in detail, the DFSZ axino interaction fll]) is stronger than the QCD 
interaction 03]) at low reheat temperature and thus the DFSZ axinos are also produced effi- 
ciently leading to a stringent mass bound or an strong upper limit on the reheat temperature 
below the electroweak scale when the axino is stable. On the other hand, if the DFSZ axino 
decays to a Higgsino and a Higgs boson, the decay occurs before thermal freeze-out of the 
usual neutralino LSP dark matter for lower values of f a . For higher values of f a , the DFSZ 
axino decays after the neutralino freeze-out and decay-produced neutralinos may overdose 
the Universe unless they re-annihilate sufficiently [l5[. For a typical choice of the param- 
eters, the neutralino annihilation rate of order 10~ 8 GeV~ 2 turns out to lead to a right 
amount of dark matter abundance, which makes a light Higgsino LSP a natural dark matter 
candidate [l6|. When the Higgsino-like dark matter contains a sizable fraction of bino, its 
spin-independent nucleonic scattering cross-section in the decoupling limit of heavy Higgs 
bosons and sfermions is shown to be of order 10 _45 cm 2 which is within the reach of future 
direct detection experiments. 



II. AXINO DARK MATTER 

A light DFSZ axino: Below the axion scale f a , the axino-Higgsino-Higgs Yukawa 
interaction (j3J) is generated and plays a major role in producing the axino number density 
from thermal regeneration. The axino number density in unit of the entropy density, Ya = 
ria/s, is determined by solving the Boltzmann equation: 

„ff§ = 7 , (5) 

where s = (27r 2 /45)(7 !| ,T 3 , H = 0.33*Jg^T 2 /Mp is the Hubble parameter, g* is the number 
of the degrees of freedom in thermal equilibrium, and z = m/T with the relevant mass 
parameter m of the process 7. For our calculation of the axino production, we will consider 
decay (or inverse decay) processes coming from the [i term interaction (j4j) as their effect is 



expected to be larger than scattering processes when decay channels are open [17 . 
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Let us first consider the decay process Hi — > aH\ taking the Higgs mass larger than 
the Higgsino mass fi before the electroweak symmetry breaking. Integrating Eq. (|5]), one 
finds 

135fl T Hl M P f°° 3 



4tt 4 £* /2 rn 2 Hl 

where g = 4 including the Higgs degrees of freedom and the final axino number, ~ 
c H(/ i //a) 2m ^i/8 7r an d z = rriHx/T. This gives 

a U l72keVy UoOGeW \ m H J\ f a /c H ) ' 10 

for 0, = 200. 

When the temperature of the Universe falls below the critical temperature, T c ~ 100 
GeV, of the electroweak phase transition, new decay channels can open. In the decoupling 
limit of heavy Higgs bosons, the coupling of the lightest Higgs boson, h, resulting from the 
interaction (]3J) is 

4 = chj- ~a[c?Hl + Sp H$\h + h.c, (8) 

la 

where the angle j3 is defined by tp = (H® ) /{H® ) . Assuming almost degenerate Higgsinos 
(with the mass « fi) decaying to a + h, one obtains the axino abundance: 



135g TfiMp 
4vrV 3/2 m \ 



^ " . ""3/2 : "J" I * Z K x {z)dz, (9) 



where z c = m^/T c , g = 2 and « c^(/i// a ) 2 m^/167r. Taking — ji — 500 GeV and 
T c = 100 GeV the integration in Eq. <M) from z c = 5 gives 



V300keW V'"/f/ V Udl I 

for gf» = 100. Another potentially important process for the DFSZ axino production comes 
from the supersymmetric top quark Yukawa coupling: 

£* = c^j- [att c + at c t] + h.c. , (11) 

J a 

where q quantifies the axino fraction in the Higgsino H2. This opens up a new channel of 
the stop decay to an axino and the resulting axino population can be calculated as in Eq. ([9]) 
with a replacement of \i — > m t , — > m^c and z c = m^ c /T c . One would get a similar 
result as in Eq. ffTUj) except a more suppression of (m ( /m t - t - c ) 2 for mjfc ~ m^. We will use 
the stop decay for the case of a heavy axino. 

The above results can be compared with the axino production from the QCD interaction 

n^ CD e « 11 (^-) (-^-] ( 10UGeV ) 2 (12) 

"a ft ~ U - i U44MevHl0 3 GeVy' V /a 7' l j 

where Tr is the reheat temperature of the Universe. This shows that the KSVZ axino pro- 
duction is more effective than the DFSZ axino production for higher reheat temperature: 
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Tfj > 10 5 ~ 6 GeV for which the small loop factor (g 2 /327r 2 ) of the QCD interaction term 
flTJ can be overcome. In any case, the axino mass must be much smaller than the elec- 
troweak scale which requires some special arrangement in the supersymmetric axion sector 
as mentioned in the Introduction. 

A heavy DFSZ axino: When the axino mass takes its typical value of the soft super- 
symmetry breaking scale in gravity mediation, the above strong bounds can be evaded if the 
axino population is diluted by entropy dumping or the axino decays to a lighter dark matter 



candidate like the gravitino [14L Il8l|. The first possibility has been discussed in the context 
of thermal inflation 13J. In this case, the axino number density is suppressed by the Boltz- 
mann factor and thus is quite sensitive to the reheat temperature after (thermal) inflation 
or any entropy dumping. To derive an appropriate reheat temperature for which the axino 
abundance saturates the required dark matter density of the Universe: finish 2 ~ 0.11, let 
us take the axino production through the stop decay ffTTl) after the electroweak symmetry 
breaking. Modifying the result OH]), the axino number density is given by 

47T 4 #* m t Jz R 

where g = 2, « c^(m t // a ) 2 m t -/87r and zr = ttii/Tr. Here the integration is essentially 
determined by the value of Zr and the integration to the infinity is justified if the stop 
freeze-out temperature is smaller than Tr. For zr 3> 1, the integration is well approximated 
by ^tt/2z 5 i ( 2 exp(-ZR). Thus, the axino relic density produced from the stop decay at low 
reheat temperature becomes 

fltfc. „ 0.1 (™) f^XV (^Lf 2 e— (14) 

V500GeV/ V rn~ t J \ f a J V20.2/ v ; 

for m t = 172 GeV, g* = 100 and c t = 1. Taking f a = 10 10 GeV, one gets zr = 25.4 implying 
that the dangerous heavy axino relic density can be suppressed sufficiently for Tr < 40 GeV. 
A stronger bound can be obtained, e.g., from the coupling (jSJ) if a heavy Higgsino decay 
channel to the axino is open. 



III. LIGHT HIGGSINO DARK MATTER 

A more interesting possibility appears when a heavy axino is allowed to decay to the 
LSP without resorting to the above option of low reheat temperature. The DFSZ axino 
can decay to a usual neutralino through the \i term coupling (jSJ) leading to overabundant 
population of the LSP dark matter. If the decay occurs before freeze-out of the dark matter, 
such a overproduction is harmless as the produced dark matter equilibrates and its relic 
density is determined by the usual freeze-out calculation Although the axino decays 
later, the overproduced dark matter population can get depleted appropriately throu gh a 
strong re-annihilation and thus can settle down to the required value of VLpMh 2 w 0.11 [15| . 

In fact, the Kim-Nilles mechanism provides a natural framework to realize such a re- 
annihilating dark matter scenario. As we will see, the annihilation rate of a light Higgsino 
LSP of order lO^GeV -2 and the DFSZ axino decay rate turn out to be in the right range 
to make a good dark matter candidate the Higgsino LSP produced from the DFSZ axino 
decay a — > hH for the axion scale of f a m 10 n ~ 12 GeV. 
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Let us first consider if the heavy DFSZ axino population can be large enough. From the 
discussions in the previous section, the axino can be produced from the stop decay as in 
( TT3|) but with the replacement of Zr — > z c « m t -/T c . This gives Y a ~ 5 x 10" 9 for g„ = 100, 
f a = 10 11 GeV, mi = ITeV and T c = 100 GeV. A more effective and inevitable process is 
the inverse decay hH — > a. From the couplings in Eq. (jBJ), we find 

135 Y~ a M P [°° 3 
4tt 4 ^ 3/2 wig 4 



1n -/ n \2 /10 n GeV\ V500GeV\ 

for = 100 and z c = m a /T c = 5. This is more than 10 4 times larger than the required 
value for a right dark matter density. If some decay mode is kinematically allowed before 
the electroweak symmetry breaking, such as a — > HH, the integration range of z — [0, oo] 
in Eq. (I15p can be taken to increase the axino population by factor of 6. Thus, we can 
safely conclude that the Higgsino dark matter population from the DFSZ axino decay is 
large enough to require re-annihilation. 

Next question is whether the decay temperature of the axino is larger than the free-out 
temperature of the dark matter. The axino decay process a — > hH (H may decay further 
to the LSP if it is not the LSP.) has the decay rate: 

2 / \ 2 



^ " Tfa [fj - (16 » 

which corresponds to the decay temperature T D = gZ^y/ST^Mp. Defining xd = rn x /Tr) 
where x denotes the LSP (the lighter Higgsino in our case), one gets 

^V/V 500GeV y /2 (m x \ ( f a /c H \ 
XD ~ M W> {-^-) {TJKT^Gev)- (17) 

Thus, requiring the decay temperature smaller than the freeze-out temperature Tf, that is, 
xd > Xf where Xf = m x /Tf w 25, one needs f a > 7 x 10 10 GeV. Having To < Tf, the 
Boltzmann equation for the LSP abundance can be solved as 15 



y-i ~ v-ifT ) 4. ^av)s(T d) 

Y x ~ Y x ™+ h{Td , , (18) 



where the first term Y x (Td) 1 w Y~ , as given in Eq. (I15p . is much smaller than the second 
term. Thus, the dark matter relic density is determined by 

W = ^gg^, (i9) 

which requires the dark matter annihilation rate: ((Ja v ) = %d x 10~ 10 GeV -2 . From the 
relation (fTTI) . we get id ~ 340 for the (approximate) upper end of f a /cn = 10 12 GeV and 
the lower value of xd is put by Xf ~ 25. Therefore, we conclude that the annihilation rate 
in the range of 

2.5 x 10~ 9 GeV~ 2 < (a A v) < 3.4 x 10~ 8 GeV- 2 (20) 
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is required for the neutralino LSP. Such a large annihilation rate can be obtained indeed 
for the Higgsino-like LSP as can be seen from the general parameter scan made in Ref. 20 
including co-annihilations. In our case, co-annihilation does not occur as the Higgsinos are 
not in thermal equilibrium and thus the heavier states decay away to the lighter ones. 

For the calculation of the neutralino LSP annihilation rate in the mass range 200 ~ 500 
GeV, we take the leading s-wave contributions in the decoupling limit of heavy Higgs bosons 



and all the sfermions. All the relevant terms can be found in Ref. [19[ and let us recollect 
them in the below. We use the following convention for the diagonalization of the chargino 
and neutralino mass matrices: 



An 



U nl W~ + U n2 H^, 
V nl W + + V n2 H+, 
N n B + N 12 W 3 + N 13 H° 



(21) 



Note that we use the notation x = Xv The s-wave contributions to the annihilation rate in 
the W /± 14 /=F , Z°Z°, Z°h and ti final states are given by 



[a A v) ww 
(cr A v)zz 
{°Av)zh 



na 



3/2 
2 HWW 

m?, 



mi 

2 «3/2 



A 

™2 Pzh 



11 



m 



x 



3/2 
ZZ 



.71=1 

" 4 

.71=1 
' 4 

E 



U 2 



V s 



+ r l± - r w 

ATI 



+ r\ - r% 



-Z n T n (m x o - 



+ 



Z m 2 



m z{l + r 2 o - r|/2 - rl/2) c w m\ 



3tt^ ffi 



1/2 



16c' 



IV 



m 



x l 



m t m x 



m. 



(22) 



where rx = mx / m x and (3xy = yl - (r x + 7~y)/2 
in the above equations are defined as follows: 



x 



/y)/4. The couplings appearing 



U n = --^=N 13 U n2 + N 12 U nl , 

V n = -^NuVra + NuVn!, 

V2 

Z n = -7V 2 3 + iV 2 4 

T n = -(c^g - s N u ){N n2 - t w N nX ) n). 



(23) 



The LSP annihilation rate is shown in Fig. 1 by blue solid contours in the \x-M\ plane 
assuming M 2 = 2M\ and tp = 10. Remark that the annihilation rates for the whole Higgsino- 
like LSP region (and the bino-like region with a large mixture of Higgsino) are in the range 
of (1 ~ 3) x 10~ 8 GeV -2 which fits well into the desired range for re-annihilation f l20p . As 
can be expected, the most bino-like LSP region with a smaller Higgsino mixture (/i 3> Mi) 
is completely excluded as the annihilation rate becomes smaller than about 2 x 10 -9 GeV~ 2 
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FIG. 1: In the \jl-M\ plane showed is the LSP annihilation rate (blue solid contours) and the spin- 
independent nucleonic scattering cross-section (red dashed contours) in the unit of 10 _8 GeV -2 and 
10 _45 cm 2 , respectively. The dotted lines show the LSP mass. The calculation is made with the 
lightest Higgs mass 115 GeV, assuming Mi = 2M\ and the decoupling limit of heavy Higgs bosons 
and all the sfermions. The nuclear form factor is taken to be fx = 0.27. 



making Q x h 2 > 0.11. In the same figure also shown is the spin-independent nucleonic 



cross-section which is calculated using the following formula [19 



&si 



m W m h 



rp2 f 2 

4 1 1JT 



(24) 



where fi n is the reduced mass, m n is the nucleon mass, and fx is the nuclear form factor 
h = fr lq +2fT G /7 with f lq = J2 g = u ,d,s h q and f Ta = 1- h lq - Note that the direct detection 
rate has a large uncertainty due to the nucleon form factor uncertainty (see, e.g., 21] ). In 
our numerical calculation, we choose fx = 0.27 adopting the values of fx u + fr d = 0.05, and 
fx s = 0.013 from the recent lattice calculation [22| . 

The nucleon-LSP scattering depends sensitively on the bino fraction in Ti and the lightest 
Higgs boson mass. For the Higgs mass rrih = 115 GeV, the nucleonic scattering cross-section 
is shown to be as large as about 2 x 10 _44 cm 2 for the Higgsino-like LSP with a large bino 
mixture in the mass range m x = (170 ~ 190) GeV. 

Before closing this section, let us make remarks on the impact of the saxion s, the scalar 
partner of the axion 23j, whose mass is expected to be order of the supersymmetry breaking 
scale ~ 10 2 " 3 GeV. In the DFSZ axion model j24[, the dominant interaction terms of the 
saxion contain 



s 

fa 



chh(H 1 H 2 + h.c.) + d H ii l 



(25) 



s 



where the Heavy Higgs bosons are decoupled. Here, x = J2i Qi v i/ fa depends on the charges 
cjj and the vacuum expectation values Vi of the Peccei-Quinn symmetry breaking fields Si 
To|. and c' H parameterizes the dependence on soft supersymmetry breaking terms. Note 



that the axion scale is given by f a = Qi v h Like the axino, the saxion population can 
be produced thermally from the second and third terms in Eq. (125]) . When the decay of 
the saxion to the Higgsinos, s — > HH, is allowed, its inverse decay determines the thermal 
saxion population as 



405 T^bhMp 
8% 3 gl /2 m l 

2 .5xW^)Y™V(M, (26) 



y. 



v200GeVy V fa/en ) \ m s ) 

with = 100 and the decay rate = c 2 H (fi/ f a ) 2 m s /ATi. If the Higgsino channel is not 

open, the saxion decay to the Higgs bosons, s — > hh, from the third term in Eq. (I2"5"j) becomes 
the main source for the thermal saxion production which leads to the saxion abundance 
comparable to Eq. (1151) . The saxion population can arise also from a coherent oscillation 
driven by an initial misalignment, sf. 



(27) 

if the reheat temperature is high enough |25_j]. When the reheat temperature is low, Tr < 10 9 
GeV, the coherent production is typically smaller than the thermal production ( 1261) . The 
contribution of Eq. 027|) can be larger than the thermal contribution if Sj ~ f a ~ 10 12 
GeV, and thus the LSP dark matter can mostly come from the saxion decay. 

Contrary to the axino case, however, the overproduction problem of the LSP dark matter 
from the saxion decay is model-dependent. First of all, there is no such a problem when 
the saxion decay to the LSP is forbidden kinematically. Although it is allowed, the saxion 
may mainly decay to the axions through the first term in Eq. (1251) giving the decay rate 
T^aa = x 2 rn? s 1 64:71 f 2 . In this case, the produced axions are red-shifted away and the LSP 
dark matter production is suppressed by the factor T s ^fj^/T s ^ aa « 16c 2 H fi 2 /x 2 m 2 s . 



IV. CONCLUSION 

Dark matter property is discussed in the Kim-Nilles mechansim which realizes the DFSZ 
axion solution to the strong CP problem and solves the \x problem at the same time. The 
DFSZ axino regeneration by the /x term interactions (which is independent of the reheat 
temperature) is more efficient than the KSVZ axino regeneration from the QCD interaction 
at the reheat temperature of order TeV, and thus puts also a stringent bound on the stable 
axino mass. The DFSZ axino dark matter mass can vary from 0(0.1 keV) and (9(10 MeV) 
depending on the allowed decay processes and the axion scale f a ~ 10 10 ~ 10 12 GeV. When 
the stable axino has the mass of order the electroweak scale, the reheat temperature has to 
be low enough to suppress the regeneration rate and, of course, dark matter can consist of 
heavy axinos if the reheat temperature is appropriately adjusted. 

The DFSZ axino may decay to a neutralino through the \i term and the overproduced 
neutralinos can be thermalized or re-annihilate sufficiently to saturate the observed dark 
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matter density. A light Higgsino-like LSP is shown to have the annihilation rate (1 ~ 
3) x 10 _8 GeV -2 and thus becomes a good dark matter candidate realizing the re-annihilation 
scenario for the axion scale f a ~ X0 llr " 12 GeV. Furthermore, the Higgsino-like LSP with a 
sizable mixture of the bino component has the spin-independent nucleonic cross-section as 
large as 2 x 10~ 44 cm 2 for the Higgs mass 115 GeV, and thus might be observed in the future 
direct detection experiments. 
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